Gapped optical excitations from gapless phases: imperfect nesting in unconventional 

density waves 

Balazs Dor^ and Andras Vanyolos 

Department of Physics, Budapest University of Technology and Economics, H-1521 Budapest, Hungary 

Kazumi Maki 

Department of Physics and Astronomy, University of Southern California, Los Angeles CA 90089-0484, USA 

Attila Virosztek 

Department of Physics, Budapest University of Technology and Economics, H-1521 Budapest, Hungary and 
Research Institute for Solid State Physics and Optics, P.O.Box 49, H-1525 Budapest, Hungary 

(Dated: February 2, 2008) 

We consider the effect of imperfect nesting in quasi-one-dimensional unconventional density waves 
in the case, when the imperfect nesting and the gap depends on the same wavevector component. 
The phase diagram is very similar to that in a conventional density wave. The density of states is 
highly asymmetric with respect to the Fermi energy. The optical conductivity at T = remains 
unchanged for small deviations from perfect nesting. For higher imperfect nesting parameter, an 
optical gap opens, and considerable amount of spectral weight is transferred to higher frequencies. 
This makes the optical response of our system very similar to that of a conventional density wave. 
Qualitatively similar results are expected in d-density waves. 

PACS numbers: 75.30.Fv, 71.45.Lr, 72.15.Eb, 72.15.Nj 

I. INTRODUCTION 

The basic ingredient of the density wave (DW) formation is a band structure consisting of a pair of Fermi sheets, 
which can be nested to each other with a certain wavevector, giving rise to the density wave instabilitjii. In real 
materials, however, this condition is not perfectly fulfilled: £(k) -I- e(k — Q) 7^ 0. In quasi-one-dimensional models 
studied during the early history of DW, one can choose it as £(k) + e(k — Q) = 2eo cos(26fcj,), which shows the 
deviation from the one dimensionalitySiS,. In conventional CDW such as NbSca, the depression of the transition 
temperature under pressure is described in terms of the pressure dependence of imperfect nesting, and the large ratio 
of 2A/Tc is also interpreted^*^. Similarly in field-induced SDW many features are successfully described by this 
modeli2^. The general consequence of eg is the destruction of the density wave phase: imperfect nesting depresses 
the DW transition temperature and destroys completely the density wave when eg becomes larger than a critical 
value. Also the imperfect nesting term gives rise to dip structures in the angular dependent magnetoresistance in 
a-(BEDT-TTF)2KHg(SCN)4'^ and Bechgaard salts (TMTSF)2PFe^. This motivates us to incorporate the effect of 
imperfect nesting in unconventional density wave (UDW) theory. UDW is a density wave, whose gap function depends 
on the wavevector, vanishes on certain points of the Fermi surface, allowing for low energy excitations. The average 
of the gap over the Fermi surface is zero, causing the lack of periodic modulation of the charge and spin density. Such 
systems have been studied and proposed over the years in a variety of systema^iifl. These include heavy fermions 
like URuzSii^i^^i", CeCoIng", organic conductors as a-(BEDT-TTF)2KHg(SCN)^ and (TMTSF)2PFf^, high 
superconductorsiSiiLi^. Two different models are possible: 2D or 3D when the gap and the imperfect nesting depends 
on the same or different wavevector component, respectively. Previously we have analyzed the properties of the 3D 
modeliSi, and now we turn to the investigation of the 2D one. 

The object of the present paper is to extend the analysis of Refs. IT^l20l to the presence of imperfect nesting when the 
gap and the imperfect nesting depends on the same wavevector component. We discuss the temperature dependence 
of the order parameter for different eo's. The phase boundary is almost the same as in a conventional DW. The 
chemical potential is shifted from its original value of the metallic state due to the presence of imperfect nesting. The 
temperature dependence of the order parameter, A(T, eg) is anomalous: although it decreases monotonically with 
increasing temperature, but exhibits a sharp cusp at A(T, eq) = 2eo- In the density of states (DOS) the particle-hole 
symmetry is broken for the 2D model, leading to asymmetric density of states with respect to the Fermi energy. For 
high values of eo, the zero of the density of states at the Fermi energy disappears, and DOS becomes finite for all 
energies. The optical conductivity is not affected by the deterioration of perfect nesting in a wide parameter range. 
By further increasing eo, the divergent peak at 2A is divided into two new peaks. Moreover, a finite optical gap shows 
up at T = in spite of the finite density of states. Similar behaviour was identified in two dimensional UDW (the 
so-called d-density wavei^): deviations from perfect nesting induce a finite optical gap^^. In clean systems, the weight 
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of the Dirac delta peak at zero frequency is finite for all temperatures. We expect similar results in d-density waves 
as well. 



II. PHASE DIAGRAM 



The single-particle electron thermal Green's function of UDW is given bjsSSiS^ 

G^^ikjujn) = iun - ^(k)/93 - picr3ReA(k) - p2(73lmA{k), (1) 

where pi and ct^ {i — 1, 2, 3) are the Pauli matrices acting on momentum and spin space, respectively, and for UCDW (73 
should be replaced by 1. A(k) — Ae*"^ cos(fefey) or sin(6fcj,). ip is the unrestricted phase (due to incommensurability) 
of the density wave, but its explicit value turns out to be irrelevant for our discussion, ^(k) is the kinetic energy 
spectrum, ujn is the fermionic Matsubara frequency. The effect of imperfect nesting is incorporated in the theory by 
replacing the Matsubara frequency in the single particle Green's function with w„ + i{eo cos{2bky) — 6iJ,)^^^, where 
(5/i is the change of the chemical potential due to the change in the spectrum. The order parametes^ is assumed 
to depend on the wavevector like A(k) = Asin(bky) or Acos(6fcj,). The second order phase boundary is given by 
eo — Ao (Tc), where Ao(T) is the temperature dependence of the gap in a perfectly nested conventional DW with 
TcQ transition temperature. Tcq is the transition temperature in the absence of imperfect nesting. This is almost 
the complete phase diagram. At high temperature when T becomes of order of eo the deviation from perfect nesting 
becomes irrelevant and the best nesting vector is Q = (2fc^, 7r/&, tt/c). In the conventional scenario two DW phases 
can occur«2Si, characterized by slightly different wave vectors and Q is replaced by a temperature dependent wave 
vector, opening a narrow region above the critical nesting at low temperatures. For the present model, the possibility 
of ordering with different wave vector is there, although its examination is beyond the scope of the present discussion. 
The critical nesting is given by eo = -\/eAoo/2 ~ 0.82Aoo, where Aqo is the gap in a perfectly nested system at 
zero temperature. The order parameter remains unchanged for eg < Aoo/2, and vanishes sharply as eo approaches 
its critical value. This together with the phase diagram is shown in Fig. ^ The most interesting consequence of 
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FIG. 1: The phase diagram (solid line) and the zero temperature order parameter (dashed line) are plotted in the presence of 
imperfect nesting. 

imperfect nesting is that the chemical potential does not remain constant under the density wave formation. Its shift 
is given by Sp = eoe(A(r, eo) - 2eo) + A(T, eo)^/(4eo)e(2eo - A(r, eo)), where e(x) is the heaviside function. This 
behaviour can readily be seen from the density of states, where for any finite A the total number of states below 
the Fermi energy is regained only by shifting the Fermi energy as given above. Note that this change belongs to a 
sinusoidal gap while for a cosinusoidal gap the sign of the shift is reversed. The change in the spectrum in the presence 
of imperfect nesting is shown in Fig. [SJ which is given by 
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In the perfectly nested case, the low energy part of the spectrum consists of Dirac cones with peaks at the Fermi 
energjiSS. For small Eq, the spectrum is still crossed by the Fermi energy at the zeros of the gap. By increasing eg, a 
broad bump develops in the upper band, and crosses the Fermi energy. At this point, a large number of possible states 
becomes available, and the chemical potential starts decreasing to keep the total number of particles unchanged. 
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FIG. 2: The evolution of the quasiparticle spectrum is shown, viewed from the direction of the kx axis, for A(k) — Asm{bky) in 
the presence of imperfect nesting for eo/A = 0, 0.3, 0.5 and 1 from left to right, top to bottom. The horizontal line denotes the 
Fermi energy. The band structure is chosen as e(k) = —2ta cos{akx) — 2tb cos{bky) + eo cos{2bky) with parameters as fa/A = 2, 
tt/A = 0.1 at half filling. 

A direct consequence of this shift is a cusp in the temperature dependence of A at A = 2eo, since at this point the 
chemical potential changes. This feature is shown in Fig. which is obtained from the numerical solution of the gap 
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equation: 



2ir 



sin^(y)dy 



n I ^ (w„ + i(eo cos(2y) - 5^)Y + sin^(y) 



(3) 



where 1/ > is the interaction responsible for the UDW formation, A'o is the density of states in the normal state at 
the Fermi energy per spin. 




FIG. 3: The temperature dependence of the order parameter for the 2D model is shown for eo/Aoo = 0, 0.2, 0.3, 0.45, 0.55 and 
0.7 from right to left. The cusp shows up only for 2eo < Aoo- 



III. DENSITY OF STATES 



The quasi-particle density of states is given by 

i ^(i? + 5/i-eocos(2y))2-A2sin2(y) 

The energy variables are expressed in units of eg, i.e. a = A/eq and x — {E + (5/i)/eo, the energy is measured from 
the new Fermi energy. The density of states is obtained as 



'^VnW p-qj \2V pq / p-q \ -^pq ^\ pq // 



for X > aV8 + l, where p = ^(m- 1)2 + ^2, q = Vm2 +^2, m = (a2 - 4(x - l))/8, n = ay'-a^ + 8{x- l)/8, K{z) 
and n(n, z) are the complete elliptic integrals of the first and third kind2^. In the remaining regions the DOS reads 
as 

N{x, a) = Noieia - 4)/i(x, a) + 9(4 - a)Mx, a)), ^ + l>x>a-l, 

o 

N{x, a) = iVoSgn(a; — l)/2(a;, a), a — 1 > a; > —a — 1, 

N{x,a) = -Nof3ix,a), -a-l>x, (6) 

where the following notations are used: 
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f2{x,a) = 



f3ix,a) = 




yi 



' iy2 - i)yi 
y2 - yi 



- 2/1)2/2 

l + 2yi+ sgn{x - l)ix + 1))K 




and 2/1 = (a2 - 4(2; - 1) - ay/a^ - 8{x - l))/8, ya = (a^ - 4(.t - 1) + ay/a^ - 8{x - l))/8. 



(8) 



(9) 



The particle-hole symmetry is broken, which can be readily seen from the behaviour of the peaks in the density of 
states, which slide from ±A to — A — eo — 6pi below the Fermi surface, while above it to A — eg — Sfi for 4eo < A and to 
£0 + A^/8eo — otherwise. Also the zero in DOS is at the new Fermi energy for eg < A/2, and for larger eg there exists 
no zero in the DOS. The density of states is plotted in Fig. 0] These statements correspond to A(k) — Asin(A:j^6), 
while for a cosinusoidal gap E —E change is needed in the density of states. 
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FIG. 4: The density of states as a function of energy is shown in the left panel for eo/A = (solid line), 0.2 (dashed hne), 0.45 
(dashed-dotted line). In the right panel eo/A = 0.55 (solid line), 0.7 (dashed line) and 0.8 (dashed-dotted hne) are used. 

The residual density of states (i.e. g2D{E — 0)) is given by ./VoO(2eo — A). Since on the Fermi surface the DOS 
vanishes in the same way for eg ^ Ago than in the perfectly nested case, the specific heat increases quadratically with 
temperature close to T = K in this region, while for large eg it equals to the specific heat in the normal state. 



IV. OPTICAL CONDUCTIVITY 



In this section we investigate the quasiparticle contribution to the optical conductivity. For simplicity we neglect the 
effect of the quasiparticle damping due to impurity scattering for example. The quasiparticle part of the conductivity 
contains relevant information about the system in the perpendicular cases {y and z) when the effect of the collective 
contributions can be neglected. The regular part of the optical conductivity (without the Dirac delta) is given by 



KeC/(-) = No4 r ^ r ^A^ne^^^^mm^ ftanh + tanh 



^(w/2)2-A2(k) V~ V 4T y ' ^ V 4r~ ' ' ' 

where WQ,(k) is the quasiparticle velocity in the a direction, , = vp, Vy = \/2tbb, = \/2tcC and rj = cq coa{2bky)—Sfi. 
From now on we restrict our investigation to the T = K case. The optical conductivity remains the same as in 
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the perfectly nested case for 2eo < A(0;eo) = Aqo- For higher eo the optical conductivity is zero for oj < G, 
G = A(\/8 — — a)/2 similarly to the effect of magnetic field where the w < 2^bH part of the conductivity is 
chopped^*, in other words a clean optical gap develops. This can readily be observed in Fig. [3 when the upper 
band crosses the Fermi energy, the chemical potential moves below the zeros of the gap, suppressing the low energy 
excitations, since only q = transitions are allowed for. Parallel to this the peak at 2A splits into 2 new peaks at 
A(-\/8 — a? + a)/2 and at A(a/2 + 2/a). For cj > 2eo(l + the optical conductivity remains unchanged compared 

to Ref. jEol The only change in the remaining region can be expressed by redefining the / functionsSS: 

2 

/(a,/3,5) = ^ [f [g^,x^ - F {g./^,x) ~ E (g^,x) + E {g^,x)) (11) 



2 

- 4q;2 



and IcosicK, (3, g) — I{a, f3, g) — Ismict, (3, g), where F{z, k) and E{z^ k) are the incomplete elliptic integrals of the first 
and second kind, x = 2/S./ujg'^ and the arguments of the / functions are obtained as 

f 1 a? auj\ 

a = max 0, - — — (13) 

V 2 8 4A y ^ ^ 

^^--^2-y + 4A'(2A) j (14) 
g = max^l,^^ (15) 

for Lj > G and for a < 2. For a>2, a — 0, f] — 1 and the / functions reduce to those in Ref. |23- Here min and max 
gives the largest and the smallest value of its arguments, respectively. With these notations the optical conductivity 
reads as 

Rea';;-'°^Lu) = e^Novl^^^^^Is^n,cos ia,P,g) , (16) 
Rea,, iu) = e'Novl^^^^^I (a, (3,g). (17) 




Lo/Aoo 

FIG. 5: The real part of the complex conductivity for the 2D model in the z direction is shown for eo/Aoo = ~ 0.5 (dotted 
line), 0.6 (dashed line), 0.7 (dashed-dotted line) and 0.8 (solid line). Note that the same curves belong to a^xi^^) by changing 
Vz to vf- The inset shows the eo dependence of the optical gap. 



FIG. 6: The real part of the complex conductivity for the 2D model for a sinusoidal gap in the y direction is shown for 
eo/Aoo = — 0.5 (dotted line), 0.6 (dashed line), 0.7 (dashed-dotted line) and 0.8 (solid line). The inset shows the eo 
dependence of the peaks. 
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FIG. 7: The real part of the complex conductivity for the 2D model for a cosinusoidal gap in the y direction is shown for 
eo/Aoo = (dotted line), 0.6 (dashed line), 0.7 (dashed-dotted line) and 0.8 (solid line). 

The optical conductivity in the three qualitatively different cases is shown in Figs. El El E| In the x direction 
the quasiparticlc part of the optical conductivity is the same as azz{'^) if we replace Vz with vp, although in the x 
direction it does not give the total conductivity since collective contributions change significantly the quasiparticlc 
part as it was shown in Ref. 29. At first sight the sum rule seems to be violated since a lot of optical weight is missing 
at small frequencies below the optical gap. But the 5{uj) part of the conductivity does not freeze out at T ^ in 
the presence of imperfect nesting, and its coefficient provides the missing area. As is well known, in the presence of 
impurity scattering, 5{lli) changes to a Drude like peak centered at w = 0. At finite temperature, the optical gap 
vanishes, but excitations below G are only possible with a probability of ^ exp(— (eg — A-^/(4eo ))/?'). 

V. CONCLUSION 

We have studied theoretically the effect of imperfect nesting in unconventional density waves. Two qualitatively 
different cases are possible: the gap and imperfect nesting depend on the same (called 2D model) or different wave 
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vector component (3D case)i^. Here we concentrated on the former. We explored the phase diagram which is 
identical to the one in conventional density wave. The zero temperature gap coefficient is not constant contrary to 
the conventional case. The chemical potential changes compared to the normal state value. The density of states 
turned out to be asymmetric with respect to the Fermi energy due to the particle-hole symmetry breaking, but the 
logarithmically divergent peaks of the eo = case remain present, but at different energies. For larger values of 
imperfect nesting (2eo > A(r, eo)), the zero at the Fermi energy disappears, and the low energy density of states 
regains its normal state form. Usually cq is thought to vary with pressure providing the opportunity to check these 
result in a wide range of parameters. The optical gap of the model in the perpendicular optical conductivity can be 
observed experimentally at low temperatures. Moreover the splitting and lowering of the resonant peak at w = 2A 
(when the wavevector dependence of the gap and the velocity coincide) or its absence (for the other kind of gap)i^ 
could provide robust signatures of the microscopic nature of the low temperature phase. 
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